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The general Lagrange-Euler formalism for the three memory circuit elements, namely, memristive, 
memcapacitive, and meminductive systems, is introduced. In addition, mutual meminductance, i.e. 
mutual inductance with a state depending on the past evolution of the system, is defined. The 
Lagrange-Euler formalism for a general circuit network, the related work-energy theorem, and the 
generalized Joule's first law are also obtained. Examples of this formalism applied to specific circuits 
are provided, and the corresponding Hamiltonian and its quantization for the case of non-dissipative 
elements are discussed. The notion of memory quanta, the quantum excitations of the memory 
degrees of freedom, is presented. Specific examples are used to show that the coupling between 
these quanta and the well-known charge quanta can lead to a splitting of degenerate levels and to 
other experimentally observable quantum effects. 



I. INTRODUCTION 

Circuit elements with memory, namely, memr 
memcapacitive and meminductive^ systems are attract- 
ing considerable attention in view of their application 
in diverse areas of science and technology, ranging from 
solid-state memoried^' to neuromorp hic ci rcuits^"!!] anc j 
understanding of biological processe d 13 * 14 !. The general 
axiomatic definition of memory elements considers any 
two fundamental circuit variables, u(t) and y(t) (i.e., cur- 
rent /, charge q, voltage V, or flux <j> = /_ V(t')dt') 
whose relation, the response g, depends also on a set, 
x = {xi}, of n state variables describing the internal 
state of the system. These v ariabl es could be, e.g., the 
spin polarization of the samp le^* 1 ^ or the position of oxy- 
gen vacancies in a thin film 1 -!]. The resulting n-th order 
it-controlled memory circuit element is described bjP 

y(t) = g(x,u,t)u(t) (1) 
x = f(x,u,t) (2) 

where / is a continuous n-dimensional vector function. 
It is assumed on physical grounds that, given an initial 
state u(t = to) at time to, Eq. (pi) admits a unique so- 
lution. If u is the current and y[t) is the voltage then 
Eqs. ([IJ, ([2]) define memory resistive (memristive) sys- 
tems. In this case g is the memristance (for memory resis- 
tance). In memory capacitive (memcapacitive) systems, 
the charge is related to the voltage so that g is the mem- 
capacitance (memory capacitance); while in memory in- 
ductive (meminductive) systems the flux is related to the 
current with g the meminductance (memory inductance) . 
These systems are characterized by a typical "pinched 
hysteretic loop" in their constitutive variables when sub- 
ject to a periodic input (with exceptions as discussed in 
Ref . H8|) . Indeed, we have recently argued that essentially 
all two-terminal electronic devices based on memory ma- 
terials and systems, when subject to time-dependent per- 
turbations, behave simply as - or as a combi natio n of - 
memristors, memcapacitors and meminductors^ 1 ^. This 
unifying description is a source of inspiration for novel 



digital and analog application d 18 | 20 | 2 jJ and allows us to 
bridge apparently different areas of researchP 

However, despite the wealth of applications and new 
ideas these concepts have generated, it is nonetheless im- 
portant to stress that so far these memory elements have 
been discussed only within their classical circuit theory 
definition, with quantum mechanics entering at best in 
the microscopic parameters t hat dete rmine the state vari- 
ables responsible for memory* 3 ! 17 l 22 l. However, it seems 
that these features are common at the nanoscale where 
the dynamical properties of electrons and ions are likely 
to depend on the histor y of the system, at least within 
certain time scales^ 3 -^. Mindful of the trend towards 
extreme miniaturization of devices of all sorts, it is thus 
natural to ask whether true quantum effects can be as- 
sociated with the memory of these systems and which 
phenomena could emerge from the quantization of mem- 
ory elements. Of course, examples of memory effects in 
quantum phenomena can be found in the specialized lit- 
erature (see, e.g., Ref. Here instead, we want to 
provide a general framework of study of the quantum ex- 
citations (memory quanta) associated to general degrees 
of freedom that lead to memory in these systems. 

We then first introduce the general Lagrange-Euler for- 
malism for these systems. This is the non-trivial exten- 
sion of the corresponding formalism for the "standard" 
circuit elements. Since it is well known that the La- 
grangian formulation of circuit elements offers great ad- 
vantages in the analysis of complex circuits^, we expect 
that this generalization would be of great value in it- 
self. Moreover, our work extends previous studies related 
to the formulation of Lagrange and Routh equations for 
non-linear circuits involving ideal memristors 2 -^ and to 
the port-Hamiltonian modeling for the case of memris- 
tive components 2 ^". In the present context our work also 
sheds light on the general relation between the internal 
degrees of freedom that lead to memory and the consti- 
tutive variables - the charge, current, voltage and flux - 
that define the different elements. Along the way we also 
define mutual meminductors, namely mutual inductors 
with memory, which add additional flexibility and hence 
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new functionalities to the field of memory elements. 

We finally proceed to quantize the corresponding equa- 
tions in the standard way. This leads us to consider 
the memory excitations of these systems. In this pa- 
per we consider only the quantization of non-dissipative 
elements, and we will devote a subsequent paper to the 
discussion of quantum effects in dissipative memory ele- 
ments. We will provide examples of applications of the 
Lagrangian formalism to selected cases and discuss ex- 
perimental conditions under which these memory quanta 
could be detected. 

This paper is organized as follows. In Sec. |TT] we in- 
troduce a general scheme of the approach. Sec. Ill is 



dedicated to the Lagrangian formulation of memristive 
systems, while Sees. |IV| and [V] deal with memcapacitive 
and meminductive systems, respectively. We then show 
how to write the Lagrangian (Sec. VI I and Hamiltonian 
(Sec. VII I of a circuit of memory elements and also give 
the work-energy theorem and generalized Joule's first law 
for such a circuit. We introduce the concept of memory 
quanta in Sec. VIII focusing on specific examples. Fi- 



nally, in Sec. |IX| we report our conclusions. 



II. LAGRANGE APPROACH 

In the Lagrange formalism, each memory circuit ele- 
ment is associated with m + 1 degrees of freedom (one 
related to a circuit variable (q or </>) and m to its inter- 
nal state (generalized coordinates, yj, j = I, ...,m)). For 
convenience, we define two multivariate vectors 



= (4>,m, -,Vm)- 



(3) 
(4) 



We note that there are two (in some cases, however, one) 
internal state variables Xi (entering Eqs. (JIJ, p])) for each 
yj. Quite generally then x = {y,y} (with y here not to 
be confused with the output variable y(t) in Eq. Q). 

A model of any particular memory circuit element con- 
sists of three components: the kinetic energy, T, the po- 
tential energy, U, and the dissipation potential W. The 
m + 1 Lagrange equations of motion are given by 



d dC 

dt dYf 



dC 



(5) 



where C = T — U is the Lagrangian, a is q or tf> 
j = 0, ...,to, and the generalized dissipation force Qy? 
is defined as 



&H 

dY a 



(6) 



While, generally, models of different memory circuit ele- 
ments involve similar terms related to internal degrees of 
freedom, the contribution from the circuit variable q or 
<p is specific for each type of memory circuit element as 
presented in the Table [TJ 



System type 


variables 


T 


u 


n 


V-controlled memristive system 


Y" 






n/M 
fly 


/-controlled memristive system 


y<t> 








V-controlled memcapacitive system 


Y" 




u$ 




(/-controlled memcapacitive system 


y<t> 


rjiC 

1 q 






(^-controlled meminductive system 


Y" 








/-controlled meminductive system 


y<t> 









TABLE I. General scheme of Lagrange description of memory 
circuit elements. Specific contributions listed in the columns 
T, U and H are given by Eqs. |f6|, @, Q, |39|, (50}, 
(571. 



The kinetic energy T may have a contribution describ- 
ing the dynamics of internal degrees of freedom and a 
specific contribution according to TableJT] The contribu- 
tion from internal degrees of freedom, T, can be written 
using symmetry arguments. First of all, since dissipative 
effects are not included in the kinetic energy, it is time- 
reversal invariant, and only even powers of yi can exist. 
In order for the transformation to canonical momenta be 
invertible, however, we must leave only quadratic terms, 
and we find T = ^.^Cijijiyj. This form, being symmet- 
ric, can be diagonalized to give 



T = T 



i2 
CiVi 



(7) 



where Cj are real positive numbers to be determined mi- 
croscopically, and is the specific contribution, if it 
exists (see Table [l]) to the kinetic energy for u-controlled 
memory circuit element, /3 = M, C or L. 

The potential energy U and dissipative potential % 
also include a specific contribution from Table [T] and con- 
tributions from internal degrees of freedom U and %: 



U = U(y,u,t) + U£, 
H = H(y,y, U> t)+Hi 



(8) 
(9) 



where y = {yi}. It is important to consider the control 
variable u as an independent parameter that can be re- 
placed by an (output) circuit variable (using, e.g., Eq. 
(fll)) only in the final equations of motion. 



III. MEMRISTIVE SYSTEMS 

There are two types of memristive systems: voltage- 
controlled and current-controlled ones^l. From Eqs. (UJ) 
and ([2]), we define voltage-controlled memristive systems 
by the equations 

I M (t) = R- 1 (x,V M ,t)V M (t), (10) 
x = f(x,V M ,t), (II) 

where Vn/(i) and /^W = <l(t) denote the voltage and 
current across the device, and R is the memristance 
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so that the Lagrangian is 



(a) 



(b) 



FIG. 1. (a) Schematic of a voltage-controlled memris- 
tive system connected to a time-dependent voltage source, 
(b) Schematic of a current-controlled memristive system con- 
nected to a time-dependent current source. 



and its inverse is the memductance (for memory conduc- 
tance). A current-controlled memristive system is such 
that the resistance a nd th e dynamics of state variables 
depend on the currenPEsl 



V M (t)=R(x,I M ,t)lM(t), 
x = f {x,I M ,t) . 



(12) 
(13) 



At this point we note that the above equations have 
been introduced to define a wide class of systems col- 
lectively called memristive^, while the name memristoi^ 
has been assigned to the ideal case of these equations, 
when R depends only on the voltage (or current) his- 
tory. Although some authors use the term memristor 
to represent any system that satisfies Eqs. ( 10 ),(|TT|) or 
(12 1,(13) we reserve this term for the ideal case onhJR 



(We will also see in Sec. VIC that such systems, like 



ideal memcapacitors and meminductors, require special 
care in the Lagrangian formulation.) We also note that, 
often, current-controlled memristive systems can be re- 
defined as voltage-controlled ones and v ice-ve rscPR In 
addition, according to Thevenin's theorerrl 29 ! 30 !, a voltage 
source V(t) in series with a resistance R is equivalent to 
a current source I{t) = V(t)/R with the same resistance 
in parallel. We could then choose to work with either 
one of these cases. However, for completeness, in the fol- 
lowing we will present the Lagrangian formalism for both 
voltage-controlled and current-controlled memristive sys- 
tems. 



A. Voltage-controlled systems 

We consider a voltage-controlled memristive system 
connected to a time-dependent voltage source V(t) as 
shown in Fig. [lja). In addition to the term due to in- 
ternal degrees of freedom discussed in Sec. [TTJ the total 
potential energy contains the usual contribution from the 
battery —qV(t), with q the charge that flows in the cir- 
cuit. There are many mechanisms for potential energy 
arising from the state variables which are affected by the 
applied bias — an example of this is the change of state 
due to electromigration (see, e.g., Ref. l24|) . The total 
potential energy is thus given by 



U = U(y,V M ,t)-qV(t) 



(14) 



C = T-U = 



E 



2 



U(y,V M ,t) + qV(t). (15) 



Here, Vm is considered as an independent parameter. 

As it is shown in Table [TJ the dissipation potential 
of voltage-controlled memristive systems includes a cir- 
cuit variable contribution Ry . We write it similarly to 
the well known Rayleigh's "dissipation potential" (for a 
constant value resistor) of the type R — Rq 2 /2, which 
gives rise to a "dissipation force" Q q = —VqR = — 
Specifically, we will use 



7-1 M 
rt-v 



R(y,V M ,t) q 2 



(16) 



At this point we stress that the memristance (as well 
as the memcapacitance and meminductance we will dis- 
cuss later) may also depend on generalized velocities, y, 
which are also included into x. This would simply modify 
the Lagrange equations of motion without changing the 
overall formalism. To simplify the notation, however, we 
will not include this dependence explicitly here, and give 
an explicit example of this case in Sec. |VB| 

For the total dissipation potential we write 



H 



R(y,V M ,t) q 2 



H(y,y,V M ,t) 



(17) 



where the last term is to be determined phenomenologi- 
cally or from a microscopic theory. 

It is straightforward to show that the equation of mo- 
tion ([|) for Y Q q = q can be written as 

V(t) = V M {t)=R(y,V M ,t)q. (18) 

This equation is of the type ( |To| . The corresponding 
equations of motion for the state variables x are 



.. , &H(y,y,V M ,t) dU(y,V M ,t) 

C ^Vi H 1 fl 

oyi <% 



= 0, 



(19) 



which show explicitly two possible physical origins of 
memristance - due to a dissipative component and/or a 
potential energy component. 

Equation (19) can be rewritten as two first-order dif- 



ferential equations of the form (11) considering both yi 
and i)i as internal state variables. Moreover, in the final 
equations we can substitute Vm by its expression in terms 
of the current q. For this purpose, Eq. ( 10 ) can be solved 
with respect to Vm- The same final procedure can also 
be used in the case of memcapacitive and meminductive 
systems considered below. 



B. Current-controlled systems 

As a simple example of a closed circuit with a current- 
controlled memristive system, we consider a source of 
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current I(t) connected to a memristive system (Fig. 
[jjb)). Here, as indicated in Eqs. (12 1 and (13), the 



output circuit variable is the voltage across the memris- 
tive system, cb — Vm (t)- That is why we use Y& set of 
variables in this case. The kinetic energy, potential en- 
ergy, and total dissipation potential in the Lagrangian 
formalism are now 



lJ2 



U = U(y,I M ,t)-<f>I(t), 

,M 



n = uf +n = 



2R(y,I M ,t] 



(20) 
(21) 

+ H(y,y,I M ,t), (22) 



where —<pl(i) is the battery term. Although 
not a necessary step (if their values are known), 
R(y,I M ,t), U(y,I M ,t), K(y,y,I M ,t) can be obtained 
from R(y,V M ,t), U(y,V M ,t) and "H(y, y, Vm, 0; corre- 
spondingly However, the solution may be multiple- 
valued in I m so that R(y,y, Im ,t) niay have multiple 
branches with the correct choice of branch depending on 
the history of the memristive system. 

The equation of motion (EOM) for <f> follows from Eq. 



for Yq = <f>, taking into account Eqs. (20), (21), and 



22|, leading to 



R(y,I M ,t)l 



M: 



which is just Eq. ( 12 ). The EOM for x is similarly found 
from Eq. ([5]), resulting in Eq. (19) except for the substi- 
tution of Vm by Im- 



C. Example 

Here we provide a specific physical example to clarify 
both the formalism and the different terms that appear 



in Eqs. (18) and (19). For this we consider a thermistor, 
namely, a temperature-de pend ent resistor. The mem- 
ristive model of thermistoipliSI utilizes a single internal 
state variable, the absolute temperature of thermistor, 
x = y = Ttherm, and can be formulated as first-order 
voltage-controlled memristive system^. Mathematically, 
the Lagrangian model of thermistor involves the following 
kinetic and potential energies and dissipation potentials: 

(24) 
(25) 

?/;': (26) 

(27) 



T 


= 0, 




U 


= 0, 




M 


_ R(y)q 2 


V 




2 ' 








H 


= y 


\c h y 



V 2 



(T env ~ y)S 



R{y) 

where R{y) — Roe^ tyl ^ y ~ 1 / T °^ is the temperature- 
dependent resistance, Rq denotes the resistance at a cer- 
tain temperature T , /? is a material-specific constant, 
Ch is the heat capacitance, 8 is the dissipation constant 



(«) 



(b) 



FIG. 2. (a) Schematic of a voltage-controlled memcapac- 
itive system connected to a time-dependent voltage source, 
(b) Schematic of a charge-controlled memcapacitive system 
connected to a time-dependent current source. 



of the thermistor 2 , and T env is the background (environ- 
ment) temperature. 

Using Eq. ^ for a circuit consisting of a thermis- 
tor connected to a voltage source V(t) (see Fig. [lja)), 
we recover the equations of the memristive model of 
thermistor^ 



C, 



I = 

dy 

di = 



iZoeWv-i/To) 



V, 



M ■ 



V 2 
" M 



(T e , 



(28) 
y)S. (29) 



Note that although other forms of potential and kinetic 
energy terms could produce the same Eqs. (28) and (29) 



(23) this particular one also satisfies the Joule's first law dis- 



cussed in Sec. VII C This puts severe constrains on the 
choice of Lagrangian. 



IV. MEMCAPACITIVE SYSTEMS 

We now consider memcapacitive systems^ (Fig. |2j), 
which — unlike memristive systems — store also energy. In 
particular, voltage-controlled memcapacitive systems are 
defined by Eqs. |l]) and ^ with u the voltage, Vc(t), 
across the memcapacitive system, and y(t) the charge, 
qc(t), stored in the device, leading to 



qc(t) = C(x,V C) t)Vc(t) ) 
x = f(x,V C) t) , 



(30) 
(31) 



where C is the memcapacitance. 

As in memristive systems, the above equations de- 
fine a large class of systems, with ideal memcapacitors 
those for which the memcapacitance depends only on 
the voltage history (or for charge-controlled memcapac- 
itive systems, only on the charge history)P. In addi- 
tion, it is often important to consider the energy added 
to/removed from a memcapacitive system, namely the 
quantity Uc = J t Vc{t)I{t)<1t which helps understand- 
ing whether a memcapacitive system is non-dissipative, 
dissipative, or active^. Of these, the non-dissipative 
and/or dissipative memcapacitive systems are the most 
interesting for potential applications, and we will there- 
fore focus here on these cases only. 
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A charge-controlled memcapacitive system is denned 
by the set of equations^ 

V c = C- 1 (x,qc,t)q c (t), (32) 
x = f(x,q c ,t). (33) 

A. Voltage-controlled systems 

The Lagrange model of voltage-controlled memcapaci- 
tive systems is based on Y q set of variables (see Table III). 
The specific contribution from the q degree of freedom to 
the potential energy is 



U 



2C(y,V c ,t) 



(34) 



Taking into account a voltage source connected to the 
system (Fig. [2|a)), the total potential energy is written 
as 



U 



2C(y,V c ,t) 



+ U(y,V c ,t)-qV(t). (35) 



Consequently, the Lagrangian is given by 



C = T-U = 



4- 2 2C(y,V c ,t) 



-U(y,V c ,t)+qV(t). 

(36) 

The dissipativc potential contains only the internal state 
variables contribution H(y, y, Vc, t). 

The Lagrange EOMs for voltage-controlled memcapac- 
itive systems have the form 



C(y,V c ,t) 



Ciyi 



= V(t) EE V C (t), 

&H{y,y,V c ,t) 



dU(y,V c ,t) 
dyi 



dy% 

V 2 dC(y,V c ,t) 
2 d Vi 



= 0, 



(37) 



(38) 



where in writing the last term in Eq. ( 38 ) we have made 
use of Eq. ([37]) . Its clear that Eqs. (|37[), (|38| are of the 



form of Eqs. (30), (31) In fact, Eq. (38 1 clearly shows 



that the memory may arise from both a conservative po- 
tential contribution as well as a dissipative one. 

Equation ( 38 ) describes an effective dynamical system 



and, together with Eq. (30), tells us that, in the presence 



of a periodic input of frequency oj, charge dynamics can 
be out of phase with the voltage across the memcapaci- 
tive system. Indeed, there might be delay in response of 
the internal state variables to the applied voltage leading 
to the above mentioned effect. Experimentally, it ca n be 
seen as a pinched hysteresis loop in the q — Vc plane^ ,s 



Here, as seen in Eqs. 



( 32 ) and ( 33 ) , the circuit variable is 

■V c (t), 



the voltage across the memcapacitive system, < 
instead of the current through it, q = Ic(t), as in voltage- 
controlled systems, c.f. Eqs. (30), (31). Consequently, 
our analysis should be based on the set (Table fil). 
The kinetic energy, potential energy, and total dissipation 
potential in the Lagrangian formalism are now 



T = T + T„ 



Ciy 2 i+T:C{y,qc,t)4>\ (39) 



U = U(y,q c ,t)-<pi(t), 
H = H(y,y,q C) t). 



(40) 
(41) 



The EOM for (j) is derived by applying Eq. ^ to Eqs. 
@, (|40|), and leading to Eq. @. The EOM for 
x is similarly obtained and results in Eq. ( |38[ ) except for 
the substitution of Vc by qc- 



C. Example 

As example of a voltage-controlled memcapacitive sys- 
tem we consider a parallel-plate capacitor with elastically 
suspended upper plate and a fixed lower plated When 
charge is added to the plate, the separation between 
plates changes as oppositely charged plates experience 
an attractive interaction. The internal degree of freedom 
of the elastic memcapacitive system is the position of the 
upper plate y measured from an equilibrium uncharged 
plate separation, d . The Lagrange model of elastic mem- 
capacitive system connected to a voltage source consists 
of the following kinetic and potential energies and dissi- 
pation potentials: 



T = 



my 

tC 



U = U^ +u 



q 2 ky 2 
2C(y) + ~Y 



qV(t), 



H = H 



jmy 



(42) 
(43) 
(44) 



supplemented by the expression for the memcapacitance, 
C(y) — Cq/ (1 + y/do). Here, m is the mass of the upper 
plate, 7 is a damping coefficient representing dissipation 
of the elastic oscillations, k is the spring constant, Co = 
eS/do is the equilibrium value of capacitance, S is the 
plate area, and e the permittivity of the medium. 

The first equation of motion is of the form of Eq. (|37l) . 



The second equation (31 1 is obtained substituting Eqs. 
(|42|)-(|44| into Eq. (l38|). Explicitly, we obtain the classi- 



cal harmonic oscillator equation including damping and 
driving terms: 



B. Charge-controlled systems 

We consider a circuit consisting of a current source and 
a current-controlled memcapacitive system (Fig. (2^b)). 



d 2 
dt 



2+1 Tt 



V 2 



Co 

2md (l + y/d 



= 0. 



(45) 



Here, u>o = y/k/m. To emphasize the similarity of Eq. 
(45) with Eq. (31) we note that Eq. (45) can be written 



6 



(a) 



(b) 



FIG. 3. (a) Schematic of flux-controlled meminductive 
system connected to a time-dependent voltage source, (b) 
Schematic of a current-controlled meminductive system con- 
nected to a time-dependent current source. 



as two first-order differential equations and the internal 
state variables are x\ — y and x-i = y. 



V. MEMINDUCTIVE SYSTEMS 

Let us finally consider meminductive systems^ (Fig. [3]) . 
A flux-controlled meminductive system satisfies the rela- 
tions^ 

I L = L- 1 {x,4> L ,t)4> L {t) (46) 
x = f(x,4> L ,t) (47) 

with L^ 1 the inverse meminductance. A current- 
controlled meminductive system is defined by the set of 
equations^ 

<f> L = L{x,I L ,t)I L {t), (48) 
x = f{x,I L ,t). (49) 

As in the case of memcapacitive systems, meminductive 
elements may represent non-dissipative, dissipative, or 
active devices. We are interested only in the first two 
types since they are the most important for technological 
applications. 

A. Flux-controlled systems 

Consider a circuit composed of a voltage source con- 
nected to a meminductive system as in Fig. [3][a). The 
circuit degree of freedom q in flux-controlled meminduc- 
tive systems is taken into account by the following con- 
tribution to the kinetic energy 



rpL _ L(y,<pL,t)g 2 
J- — ^ 



(50) 



The contribution to T, U and % from internal state de- 
grees of freedom are written in the general form (Eqs. 
'©'©)• Consequently, taking also a voltage source 
in Fig. pMa) into account, the Lagrangian and dissipative 
potential are written as 



2 



■L(y, <i>L,t)- 



H = H (y,y,<f>L,t) 



U(y,<j> L ,t)+qV(t), (51) 
(52) 



The EOM for the q degree of freedom is 

d(L(y,<j) L ,t)q) 



dt 



V{t) = V L (t). 



(53) 



Integrating this equation in time assuming that <p(t 
—00) = we find 



L(y,<t>L,t)q 



dt'V L {t') = ct> L {t), (54) 



which is Eq. ( 46 1 



The EOMs for the state variables are written as 

dH (y,y,<f> L ,t) 



CiVi 



diji 

dU(y,<t> L ,t) dL{y,<f> L ,t) 



2i 2 % 



0, (55) 



which again, since x = {y, y}, can be written in the form 
of Eq. (|47l 



B. Current-controlled systems 

When one considers circuits involving current- 
controlled meminductive systems (and current sources 
instead of voltage sources) , one should employ the set 
of variables (see Table [l| . Let us then consider a simple 
circuit composed of a current source directly connected 
to a meminductive system (Fig. [3|h)). The contribu- 
tion from the circuit degree of freedom <f> comes from the 
potential energy Uj term. The kinetic energy, potential 
energy, and total dissipation potential in the Lagrangian 
formalism are now 



T = T= l^Ciiji, 

i 

U = U + U? -4>I = U(y,I L ,t) + 
H=H(y,y,I L ,t). 



2L(y,I L ,t) 



(56) 
07(t), (57) 



(58) 



The EOM for 4> is the same as Eq. (|48j. The EOMs for 
yi are similarly derived and result in Eq. ( 55 ) except for 
the substitution of 4>l by II- 



C. Example 

We here provide an instructive example of an effective 
meminductive system consisting of an LCR contour in- 
ductively coupled to an inductor (Fig. [4]). In this scheme, 
the two inductors L\ and interact with each other 
magnetically. From the point of view of the voltage 
source V(t), the total system can be seen as a second- 
order flux-controlled meminductive system described by 
the general equations (46 H 47 1. The charge on the ca- 



pacitor C and the current through the inductor L2 play 
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FIG. 5. Symbol for a mutual meminductor. 



FIG. 4. Flux-controlled meminductive system based on the 
inductive coupling of a coil L\ with a LCR contour. Here, the 
mutual inductance M is equal to k\JL\L<2, where < k < 1 
is the coupling coefficient. 



the role of internal state variables. It is convenient to 
select y — qc- Consequently, Il 2 = y. 

We start by considering the circuit presented in Fig. [4] 
using the Lagrange formalism for usual circuit elements. 
The circuit is described by: 



T=^L 1 q 2 + ^L 2 y 2 + Myq, 



U 

T~Lq 



2C ' 
0, 

2 



qV(t), 



The EOMs for x and q are then found to be 
L\q + My - V(t) = 0, 



Li 



M 2 



M 



V(t) 



c 



Ry = 0. 



(59) 

(60) 
(61) 

(62) 



(63) 



(64) 



One can easily verify that Eqs. (|63|) and (|64|) describe 
the electric circuit from Fig. |4j 



Next, integrating Eq. ( |63| (the constant of integration 
is taken to be zero), we can rewrite it in the form 

which shows that the meminductance L depends on the 
generalized velocity x = y. 



D. Mutual meminductance 

After the generalization of self-inductance to memin- 
ductance, one wonders if mutual-inductance can be gen- 
eralized to memory situations as well. We consider two 



coupled inductors as in Fig. g but now assume the mu- 
tual inductance to have memory. Here, we want to de- 
scribe that part of the memory that cannot be included 
in two (self-)meminductive systems. This memory can be 
stored in the medium between the inductors with a state 
affected by the two magnetic fluxes of the inductors. It 
could also be stored in the geometry of the system by 
having, e.g., two elastic coils that can either attract or 
repel each other. Since the memory mechanism does not 
belong solely to one inductor, the relation M = k\fL\L~2, 
applicable to mutual inductance of two coils, does not ap- 
ply for mutual meminductance: k is not generally a con- 
stant independent of Lx, x and possibly some other 
parameters. 

In analogy with Eqs. (46) and (47 1, we then define a 



flux-controlled mutual meminductive system via the fol- 
lowing set of equations: 



qi = M~ 1 (x,4>Ml,4>M2,t)4>M2, 
q 2 = M~ 1 {x,4>Ml,4>M2,t)4>Ml, 
X = f(x.(l) M1 ,4> M 2,t), 



(66) 
(67) 
(68) 



where M(x,(j)Mi,4'M2,t) is the mutual meminduc- 
tance, 4>mi is the magnetic flux defined by 4>m\ = 
/_ VMi(t')dt' , VMi{t) is the voltage on the first induc- 
tor (4>M2 and V2u(t) are similarly defined), and <ji and q 2 
are the currents in the first and second inductors respec- 
tively. The circuit symbol we propose for this element is 
shown in Fig. [5] 

Regarding the Lagrangian formulation, the additions 
to the kinetic energy, potential energy and dissipation 
potential as a result of introducing this memory element 
are 

T = ^ ^2 Ci^i + M( x > 4>Mi,4>M2,t)q 1 q 2 , (69) 

i 

U = U(x,<f>MX,4>M2,t), (70) 
H = 'Hx{x ) X,<j)MX,<f>M2,t). (71) 

The corresponding current-controlled mutual memin- 
ductive systems are instead defined by the set of equa- 
tions 

<Pmi = M(x, Imi,Im2, t)I M 2, (72) 
4>M2 = M(x, hn,lM2, t)I M i, (73) 
x = f{x, Imi,Im2, t), (74) 
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where lui is the current in the i-th inductor, and 
M(x, Imi-, Im2-, t) is the mutual inductance that can be 
obtained by plugging in (j) M1 (I M i) and (j) M2 {lMi) in 
M{x,4>Mi,4'M2,i)- The Lagrangian formulation of this 
system is given by 



(75) 



U = U(x, MI M2 , MI M1 ,t) + „,t M t 2 it , (76) 

n = n x (x,x, MI M2 , M I M1 , t) , (77) 

where both U and are the same functions as defined 
in the above flux-controlled case. 



VI. LAGRANGIAN OF A GENERAL CIRCUIT 

We now have all the ingredients to write down the 
Lagrangian for a general circuit network composed of an 
arbitrary combination of memristive, memcapacitive and 
meminductive systems and their standard counterparts. 
These circuits may be powered by an arbitrary set of 
voltage sources Vk(t) (Sec. VIA), for which the fluxes 



4>k (t) are defined as, e.g. , in Eq. ( 54 1 , or by a set of current 
sources It(t) (Sec. VI B ). When both voltage sources and 
current sources are present, one can convert the latter to 
the former using the Thevenin's theorem, or the former 
to the latter using the Norton's theorerrPH, thus ensuring 
only one type of power source is present. Below we briefly 
outline the recipe to write the Lagrangian of a general 
circuit for both cases. 



A. Circuits with voltage sources 

A general electronic circuit powered by voltage sources 
can be described as a combination of I indivisible loops, 
i.e., ones that do not contain internal loops. Within the 
j-th (j — 1, . . . , I) loop one should consider the charge 
qj as the circuit variable and take into account general- 
ized coordinates of elements involved in this loop. For 
simplicity, we rename the generalized coordinates for the 
whole circuit as yi (i = 1, . . . , k). 

The current in each branch of the circuit is the sum of 
contributions from indivisible loops it belongs to. Using 
this fact, we can write the Lagrangian for each element 
in the branch. The element's Lagrangian is taken in the 
voltage-controlled form for memristive and memcapaci- 
tive systems and in the flux-controlled form for memin- 
ductive ones. 

The sum of the Lagrangians of individual elements of 
the circuit gives the circuit's Lagrangian, while the sum 
of the dissipation potentials gives the circuit's dissipa- 
tion potential. The circuit's Lagrangian and dissipation 
potential depend on qj , qj , yi and yi and result in I + k 
EOMs. The EOM obtained for qj gives Kirchhoff 's volt- 
age law (KVL) for the j-th loop because of the linear- 



ity of the Euler-Lagrange equations, and because each 
component in the loop was shown above to give the cor- 
rect voltage term. Kirchhoff 's current law (KCL), on the 
other hand, is automatically satisfied by the choice of 
loop current variables. 



B. Circuits with current sources 

When a circuit is powered by current sources, the cir- 
cuit variable is the flux <fij in the j-th (j = 1, . . . , I) junc- 
tion, while 4>j is the electric potential at the junction. 
Using this definition, the flux or voltage across each ele- 
ment in the network can be found via the difference of the 
fluxes or potentials in the junctions at its ends, enabling 
one to write the element's Lagrangian and dissipation po- 
tential in the current-controlled formalism for memristive 
and meminductive systems and in the charge-controlled 
formalism for memcapacitive ones. As in the voltage- 
controlled case, the circuit's Lagrangian or dissipation 
potential is the sum of the circuit element's Lagrangians 
or dissipation potentials, respectively. 

If we denote again the internal degrees of freedom of 
the whole circuit as yi (i = 1, . . . , k), we have a circuit's 
Lagrangian and dissipation potential that depend on 4>j , 
yi and fa and result in I + k EOMs. The EOM ob- 
tained for 4>j gives KCL for the j-th junction due to the 
linearity of the Euler-Lagrange equations, and because 
each element ending on the junction was shown above to 
give the correct current term. KVL, on the other hand, is 
automatically satisfied by the choice of junction potential 
variables. 

This formalism has a complementary nature and can 
be viewed as the dual formalism to that for circuits with 
voltage sources. This conclusion will be reinforced in Sec. 
|VII[ where we will show that the canonically conjugate 
momenta of the voltage-controlled and current-controlled 
formalisms to be fluxes and charges, respectively. 



C. Lagrangian multipliers 

When a circuit, voltage-controlled or current- 
controlled, has additional constraints — missing from the 
EOMs — relating state variables to circuit variables, the 
form of their dependence should be added to the La- 
grangian. This is achieved by the method of Lagrange 
multipliers. In particular, the Lagrange multipliers are 
convenient for describing ideal memory circuit elements 
such as the ideal memristor 18 , in which the state vari- 
able y equals the charge q flowing through the device. 
Examples for such constraints appear in Subsec. |VI D| 

D. Examples 

Consider two inductively coupled RC circuits as shown 
in Fig. [6] Initially, both capacitors are not charged and 
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FIG. 6. Schematic of flux-controlled inductively-coupled 
charging circuits with a DC voltage source 



there are no currents in the circuits. The circuits are 
coupled via mutual inductance M, which results in peri- 
odic charging and discharging of the right side capacitor 
as will be seen below. We take a memristor with R(y) 
(y = q) qualitatively similar to the one fitted recently to 
experiments on Ti02 thin films^, namely of the form 



R(y) = Ro 



R 



off 



Ro 



l + V 2 /q 2 o 



(78) 



where i? n, -Roff and go ar e parameters defined for each 
memristor. The resistance is seen to decrease from R ff 
to R on as charge flows through the memristor. Denoting 
the left and right loop charges as q\ and q 2 , respectively, 
and applying the results of Sections III IV and[Vj we ob- 
tain the following Lagrangian and dissipation potentials 
for the network, 



C 



L 



Hy - gi) 5 



ml 



2C 



(79) 
(80) 
(81) 



where A is the Lagrange multiplier corresponding to the 
circuit holonomic constraint, y — q\ (a different constrain 
would lead to a different corresponding term in Eq. (79l). 



These expressions lead to the following EOMs for the 
total system 



Lq\ - Mq 2 + R{q\)qi 
Lq 2 



C 



-v 1 = o, 



Mq\ +Rq 2 + ^ = 0, 



(82) 
(83) 



where the EOM for y, giving A = 0, and the EOM for 
A, giving y — qi, were substituted. The solution of these 
EOMs for certain values of the parameters is shown in 
Fig. [7j We note that the insertion of the memristor pro- 
duces an almost constant current instead of an exponen- 
tially decreasing one in the left loop of the circuit. The 



stabilization of the current is achieved by the decline in 
the characteristic charging time R(qi)C as the capacitor 
is charged. The memristor also modifies the exchange of 



8.xl(T 7 




-2.x icr 13 



2.X10" 7 4.X10" 7 6.X10" 7 8.X10" 7 
t (sec) 



FIG. 7. (Color online) Graphs of gi(f) (top) and £72 (i) (bot- 
tom) for the system described in Fig. p] The solid line cor- 
responds to the behavior with the memristor and the dashed 
line to the behavior with the memristor replaced by a nor- 
mal resistor of resistance R s- The parameters used were 
Vi = IV, R = 10kn, Ros = lOOkfi, R on = WOQ, C = 3pF, 
M = L = 0.3mH, qo = 10 _12 C. The initial conditions were 
set to no charge or current in any of the circuit elements. 



energy between the two circuits, giving pronounced os- 
cillations in the charge of the right loop of the circuit, 
which are absent when the memristor is substituted with 
a normal resistor. 



As a second example, consider the circuit shown in Fig. 
[8] The circuit consists of a 3 x 3 network of memristors 
connected to a voltage source. Each memristor has resis- 
tance R(yi), where j/i is the cumulative charge that flows 
through the memristor. The indivisible loops charges are 
denoted by (ft. The EOMs can be readily obtained from 
the Lagrangian and dissipation potential of the circuit 
that read 



C = qV(t) - Ai(t/i - qy) - A 2 (y 2 - 92) - A 3 (y 3 - q + q{) - X 4 (y 4 - qi + (72) - A 5 (y 5 - q 2 ) 
-A 6 (y 6 - 93 + q%) - A 7 (y 7 - q 4 + q 2 ) - A 8 (t/ 8 - q + q 3 ) - X 9 (y 9 - q 3 + q 4 ) - \w{yw 



94) 
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R(y 2 ) 



Si 



«(v„) 



FIG. 8. Schematic of a voltage-controlled 3x3 network of 
memristive systems with a voltage source. 



-Aii(z/n -qo + q 3 ) - Ai 2 (yi2 - qo + q±), (84) 
H = \{R(Vi)ql + R{V2)ql + R(y 3 )(% - qif + R{Vi)(qi - qif + R(ys)ql + i?(y 6 )(<73 - qif + i?M(<Z4 - (ft) 2 + 
R{y&)(qo - q 3 ) 2 + i?(y 9 )(<Z3 - qif + R{yw)ql + R(yn)(qo ~ qzf + R(Vi2)(q - qif}-, (85) 



where the Aj's are the Lagrange multipliers. These two 
functions can be easily generalized for the case of a N x M 
network of different memristors, greatly facilitating the 
attainment of the EOMs, the solution of which can be 
used to solve, e.g., optimization problems such as mazes 
in a massively parallel way^S. 



finite network of capacitors and inductors (see, e.g.J^j). 
Since the discussion of dissipation in Hamiltonian dy- 
namics would require an extensive treatment by itself, 
here we limit our analysis to non-dissipative systems, and 
(except for the work-energy theorem discussed below) 
leave the Hamiltonian formalism of dissipative memory 
elements for a future publication. 



VII. HAMILTON FORMALISM 

The counterpart of the Lagrange formalism is the 
Hamilton one, which is also generally the starting point 
for quantization. For non-dissipative systems, one can 
easily transform the Lagrangian to the Hamiltonian. In 
the presence of dissipation instead, this task requires par- 
ticular care. 

Dissipation is the result of the tracing out of certain 
degrees of freedom resulting in an effective (reduced) de- 
scription of the system of interest in interaction with 
these degrees of freedom. However, the microscopic pro- 
cedure of tracing out these degrees of freedom is most 
of the time difficult to carry out exactly, and dissipa- 
tion is then introduced with physically plausible " ad hoc" 
strategies. 

There are several ways to add dissipation at the level 
of circuit Hamiltonians which range from complex La- 
grangians (resulting in complex Hamiltonians p21 to the 
addition of linear dissipative elements modeled by an in- 



A. Canonically conjugate momenta 

Consider a non-dissipative network of memory ele- 
ments. In order to write the Hamiltonian, we need to 
determine the momenta pj canonically conjugate to the 
variables qj. These momenta are defined for circuits of 
voltage-controlled elements by 



Pj 



dC 



(86) 



with the same definition for circuits of current-controlled 
elements except for qj being replaced by tfij. Looking 
at the expressions for the Lagrangians of the memory 
elements discussed above, one easily finds the physical 
meaning of pj. In voltage-controlled circuits pj is the 
total of the fluxes generated by the inductors in the j- 
th loop, while in current-controlled circuits, it is the 
charge in the j-th junction. In addition, if we define 
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the canonically conjugate momentum to the internal de- 
gree of freedom yi as Zj, we readily find that for both 
voltage-controlled and current-controlled circuits 



dm 



(87) 



The Hamilton's equations for voltage-controlled cir- 
cuits then read 

* = (88) 

» = (89) 



with <7j replaced by ^ for current-controlled circuits. 

Using the results for the canonically conjugate mo- 
menta, we see that for voltage-controlled circuits Eq. 
(88) gives the current in the j-th. loop in terms of the 
magnetic flux in the inductors in each loop, while Eq. 



( 89 1 gives the change in the magnetic flux in the inductors 



in the j-th loop in terms of the charges in the loops. For 
current-controlled circuits, on the other hand, Eq. (88) 



gives the potential in the j-th junction in terms of the 
charges in the circuit junctions, and Eq. (89) gives the 



current flowing into this junction in terms of the fluxes in 
the circuit junctions. The EOMs obtained here - while 
representing the same physics - are distinctly different 
from the ones in the Lagrangian formalism and therein 
lies their value. 

With these results in mind, the Hamiltonian is de- 
fined as the Legendre transformation of the Lagrangian, 
namely 



(90) 



for voltage-controlled circuits, and with qj substituted 



by 



for current-controlled circuits. Since the kinetic 



energy in both cases is quadratic in <jj (<f>j for current- 



controlled circuits) , it is easy to see that Eq. ( 90 1 reduces 
to H = T + U. 



B. Work-energy theorem 

If we consider an arbitrary circuit with a number of 
voltage sources Vkj (of the fc-th voltage source in the j- 
th loop), we can define the work done by these sources, 
at any given time in an interval of time dt, on infinites- 
imal charges dqj in each of the indivisible loops. The 
total work done by all sources (which is not an exact 
differential) is then 



6W = J2 v kjdqj 



(91) 



For non-dissipative circuits all this work goes into the 
variation of the internal energy dE which can be com- 
puted from T + U by subtracting the contribution from 



the voltage sources, 
case thus reads 



The work-energy theorem in this 



SW = dE . 



(92) 



In the presence of current sources the work done is 



SW = J2l kj (d(f> k 



(93) 



where I k j is the current of the source between the k-th 
and j-th junctions (or if none such source exists), and 
d(f>k — d<j)j is the difference in flux on the two sides of the 
source. This work accounts for the change dE of inter- 
nal energy which derives from T + U by subtracting the 
contribution from the current sources to give a balance 
formally equal to Eq. ( 92 ) . 



C. Generalized Joule's first law 

In the dissipative case on the other hand we need to 
take into account that part of the work done by the volt- 
age sources that goes into a "generalized heat" which 
accounts for the heat generated in the resistances (if 
present) and the "heat" generated from the dissipative 
components of the state variables. 

Mathematically, this amounts to 



6W-dE = J2 ^fV dt + J2 , .'/-' 



(94) 



for voltage-controlled circuits. On the other hand, in the 
presence of current sources, Eq. (94) needs to be changed 
into 



SW- dE 



dt 



(95) 



Eqs. (94) and (95) are the generalized Joule's first laws 



for dissipative systems and are important yardsticks, to- 
gether with the Euler-Lagrange EOMs, to test the valid- 
ity of a given Lagrangian formulation. We note in partic- 
ular, the identification of the energy loss due to memory, 
given by the last terms on the right-hand side of Eqs. 
( 94 ) and ( 95 1 , which are not present in the formulation 



of standard circuit elements. 



VIII. QUANTIZATION 

Having shown the Hamiltonian formulation for classi- 
cal non-dissipative circuits, we now embark on the quan- 
tization of these Hamiltonians, which will be of impor- 
tance at low temperatures and mesoscopic/nanoscopic 
length scales. Instead of proceeding with a general cir- 
cuit, in this case we find it more instructive to first work 
out explicit examples. We consider first a voltage source 
connected in series with a memcapacitive system and a 
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meminductive system. Then we look at a current source 
connected in parallel with these systems. These two cir- 
cuits can be realized experimentally (see e.g., Ref. l3"S"j) 
and are therefore ideal test-beds for the concept of mem- 
ory quanta, namely quantized excitations of the memory 
degrees of freedom of these circuits. 



A. Example: series LC circuit 

We consider a voltage source connected in series with 
a memcapacitive system and a meminductive one as de- 
picted in Fig. [9](a). The case of such a c ircuit with no 
memory was quantized in previous work:p£12£J anc [ w \\\ 
be generalized here. The Hamiltonian for this circuit is 
found using Eqs. ^ and (35) for the memcapacitive sys- 
tem and Eqs. ( p)0| ) and (J8J) for the meminductive one, and 
reads 



-1~2 



c 2i Z 2i + 



2C{ Vl ,V u t) 2L(y 2) <M) 



+ E/i(yi,Vi,t) + U 2 {y 2 ,4>,t) - qV(t), (96) 



r 



where the index 1 corresponds to the capacitor and the 
index 2 to the inductor, y^ is the i-th memory coordinate 
of the j-th memory element, and Zji is its canonically 
conjugate momentum as defined in Eq. (87 1, q is the 
charge flowing through the circuit, and <f> is the flux on 
the inductor. V\ is the voltage on the capacitor and V(t) 
is the voltage of the source. 

Under reasonable assumption of stability of the values 
of 2/1 and y 2 , we expand Ui(yt,Vi,t) and ^2(2/2, 02, t) 
at their minima with respect to j/i and y 2 , respectively. 
Several such minima can exist for y\ or y 2 in certain 
memory elements^. In this case we should choose one 
minimum based on the initial conditions. The definitions 
of t/j are shifted by constants to make them zero at their 
respective minima. This shift does not affect the form 
of the other terms in the Hamiltonian. In addition, we 
define A(C -1 )(yi, V\,t) and A(^ 1 )(y 2 , </>, t) by 



A(C~ 1 )(yi,Vi,t) = 



1 



1 



(97) 



c{ Vl ,v x ,t) c (yi,ty 

A^- 1 )^, 0, t) = 1 - — (98) 
L{y 2 ,(p,t) L ((p,t) 

where for brevity we have defined Co(Vi, t) = C(0, Vi,t) 
and Lo(p,t) = L(0,p,t). 

Discarding constant terms in Ui and neglecting higher 
order terms in Xi, we can write Ui rs £V and 
the Hamiltonian takes the form 



H = 


H q + H x + Hint , 


(99) 


H q = 


4> 2 + q 2 qV(t) 
2Lq 2Cq 


(100) 


Hy = 


2 c 3i z j* + 2 ^ ^ J '^*' 

ij ij 


(101) 


Hint = 


AiL- 1 )^ A(C- 1 )g 2 
00' 


(102) 



where the Hamiltonian was divided into a "charge" part, 
H q , "memory" part, H y , and the "interaction" part, 

Hint- 



We next introduce the bosonic creation and annihila- 
tion operators defined by 



Lqujo 

2h 



(<?• 



LqLOq 



) 



J 3 1 



3* 



2h Lqojq 

CjiLOji . IZji 

~r)*r~ ivji + — — 
2n CjiOjji 



CjiLOji 

2h 



(103) 
(104) 
(105) 
(106) 



where and a (6^ and bji) create and destroy charge 
(memory) quanta, respectively. 

The frequency of the charge oscillator, u>q, is the circuit 
resonance frequency, (LoCo) -1 ^ 2 , while the frequencies of 
the memory quanta oscillators are analogously given by 



(107) 



Plugging these relations into the Hamiltonian in Eq. ( 99 ) 
finally gives the quantized form 



H q = huj (a^a + -) - 



2Lq(jJq 



(a + <J)V(t),(l08) 



H x = ^2%uj 0l {b) i b ji + i), 



(109) 



Hin 



UL uj 



A(L- 1 )(y 2 ,0,t)(a-at)2 



4L w 



A(C- 1 )(vi,Vi,t)(a + a i ) 2 , (HO) 



where V\ is quantized by solving the equation V\ — 
q/C(yx,Vi,t) to give V\ = V%(yi,q,t) which translates 
to V\ — V\ (yi , q, t) after quantization. 

It is now clearly seen that H q includes only terms cor- 
responding to the charge quanta, while H y includes those 
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(b) 



FIG. 9. (a) Schematic of a series voltage-controlled mem- 
ory element LC circuit, (b) Schematic of a parallel current- 
controlled memory element LC circuit. 



corresponding to memory. Hi nt couples the two quanta 
with a coupling term that has at least three ladder oper- 
ators. 

A simple example^ that illustrates this result is a cir- 
cuit involving a normal inductor of inductance L in series 
with a memcapacitor that has its upper plate of mass m 
hanging on a spring with spring constant k (Fig 
This memcapacitor could be a representation of, e.g 



10) 



nano-electromechanical systerrPHHini if the displacement 
of the upper plate from its equilibrium position is de- 
noted by y and its distance from the lower plate at this 
position is do, the capacitance can be easily seen to be 
given bjff^l 



C{x) 



C 



1 + y/d ' 



(ill) 



where Co is the capacitance at equilibrium. Using the 
above formalism to quantize the Hamiltonian and keep- 
ing only energy conserving terms, we find Eq. (110) is 
reduced to 



H,, 



4V2 



d^iLCo^^imk^^ia^b + tfa 2 ). (112) 



This type of interaction is of the same kind as the one en- 
countered in the quantum treatment of second harmonic 
generation in optics^. 

If the circuit can be built to satisfy (fc/m) 1 / 2 = 
2(LCo) -1 / 2 the interaction will produce a splitting of the 
degeneracy of the levels that is of first order in Hi nt and 
which may be large enough to be detected experimen- 
tally. (See also the Conclusions for an order of magnitude 
estimate of when to expect quantum effects to dominate.) 



B. Example: parallel LC circuit 

As a second example let us consider a parallel memory 
LC circuit as plotted in Fig. [9jb) . In this circuit a cur- 
rent source is connected in parallel with a memcapacitive 
system and a meminductive system. The Hamiltonian for 
this system can be found utilizing Eqs. (39) and ( |40| for 
the former, and Eqs. ( 56 ) and ( 57 1 for the latter resulting 



H = T + U 



1 

2 ^ 



1 

2 ^ 



c 2i z 2i ' 



2C( Vl ,4>,t) 2L(y 2 ,I 2 ,t) 



+ U 1 (y ll C- 1 ^t) + U 2 {y 2 ,LI 2 ,t)-(bI(t), (113) 



with the same definitions as in Eq. (96), except for 



being the flux in the inductor and q being the charge on 
the capacitor. I 2 is the current through the inductor and 
I(t) is the current of the source. 

Proceeding in a completely analogous way to the pre- 
vious subsection with the definitions of the ladder oper- 
ators in Eqs. (103)-(106) modified by the substitutions 
L a — > Co and q — >• (f>, we find the quantized Hamiltonian 
of this system to be 



H — H q + Hy + H int 



(114) 



H q = huo{a)a + \) - \j ^^{a + <J)I(t), (115) 



h 



(116) 



Hint 



hC u>o 



A(L- 1 )(S 2 ,/ 2 ,i)(a- a t)2 



4C u; 



A(C- 1 )(J 1 ,g,i)(a + a t)2 j (117 ) 



q(0n 
-q(t) c 



FIG. 10. Elastic memcapacitive system connected to a nor- 
mal inductor. 



where I 2 is quantized by solving the equation I 2 — 
(j)/L{y 2 ,I 2 ,t) to give I 2 = I 2 (y 2 ,cj),t) which reduces to 

I2 = -^2(2/2,0,0 after quantization. This Hamiltonian is 
very similar to the one obtained for the series LC cir- 
cuit. We will now show these two Hamiltonians to be 
the basic building blocks for the quantized Hamiltonian 
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of a general circuit with non-dissipative elements. 



C. General circuit 

We now proceed to find the quantized Hamiltonian for 
a general circuit network of memcapacitive systems, me- 
minductive systems and voltage sources. Such a circuit 
can be divided into indivisible loops each with charge qk 
as noted in Sec. |VI| Using the methods of that section 
to find T and U, one can write the Hamiltonian for the 
network as H = T + U, which, after the substitution of 
qk and 4>k with ladder operators using Eqs. ( 103 )-( [T06 1, 
reduces, apart from the interaction part, to a bilinear 
combination of them which is known to be exactly diag- 
onalizable by, e.g., a linear canonical transformation. 

With regards to the quantization of the Hamiltonian 
of current-controlled circuits, the process is similar. We 
denote the flux in each junction of the network with fa 
with a corresponding qk being the charge in the junction. 
Using the methods of Sec. |VI[ we write the Hamiltonian 
H = T + V and then quantize it by writing fa and Pk in 
terms of ladder operators using the transformation from 
the previous subsection. Like in the voltage-controlled 
case, the non-interacting part of the Hamiltonian is again 
bilinear and can be diagonalized. 



IX. CONCLUSIONS 

To summarize, in this work we introduced the general 
Lagrangian formulation for the three basic memory el- 
ements: memristive, memcapacitive and meminductive 
systems and defined a fourth memory element, a mu- 
tual meminductive system, for which we also gave the 
Lagrange formalism. We showed how to write the La- 
grangian for a general circuit, including one with current 
sources. The examples given for the Lagrangian formal- 
ism demonstrated that writing the Lagrangian and dissi- 
pation potential should be the preferred choice for finding 
the EOMs of large memory element networks. 

The Hamiltonian formalism for electric circuits was 
also generalized to include memory, although only for 
non-dissipative elements. As in previous worka^fil, we 
have found that the canonically conjugate momentum 



of charge is the flux and vice versa. The Generalized 
Joule's first law was given for general circuits including 
ones with memory elements. This law can be used to 
verify the correctness of a given Lagrangian formulation. 
Lastly, we presented a scheme for the quantization of a 
general non-dissipative memory element circuit. 

The quantum treatment of memory elements, and in 
particular the example given in the text of a memcapaci- 
tor in series with an inductor (Fig. 10 ), begs the question 



of under which conditions one can measure quantum ef- 
fects in these systems. For quantum effects to be easily 
measurable, both the thermal fluctuation energy and the 
width of the energy levels should be smaller than the os- 
cillator energy quantum^, i.e., fc^T -C hui and Q> 1, 
where Q — lu R/L is the quality factor of the oscillator, 
R the loop resistance, and L is the inductance. Possible 
values for the capacitance and inductance in mesoscopic 
circuits can be taken to be 10~ 15 F^and 10" 10 tP, re- 
spectively. If one assumes a temperature of T = 20 mK 
and circuit resistance of 10 or less, both conditions 
mentioned above are satisfied. As noted for the example 
above, the degeneracy condition, satisfiable by a memca- 
pacitoi^, will lead to an experimentally detectable split- 
ting of the degenerate energy levels as a result of the in- 
teraction between the memory quanta and charge quanta. 

Future research in this field may include extending the 
Hamiltonian formalism to dissipative circuits. One way 
to do this is, e.g., via a path-integral formulatiorP2 of 
memory elements. Along a parallel line, we expect the 
Lagrangian formalism discussed here to be of great value 
in the analysis of complex networks with memory, which 
offer both fundamental and applied research opportuni- 
ties. 
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